CYCLIC TO RANDOM TRANSPOSITION SHUFFLES 

ROSS G. PINSKY 

Abstract. Consider a permutation ct € as a deck of cards numbered 
from 1 to n and laid out in a row, where Gj denotes the number of the 
card that is in the j-th position from the left. We define two cyclic 
to random transposition shuffles. The first one works as follows: for 
j = 1, • • ■ ,n, on the j-th step transpose the card that was originally 
the j-th from the left with a random card (possibly itself). The second 
shuffle works as follows: on the j-th step, transpose the card that is 
currently in the j-th position from the left with a random card (possibly 
itself). For these shuffles, for each b G [0, 1], we calculate explicitly the 
limiting rescaled density function of a;,0 < x < 1, for the probability 
that a card with a number around bn ends up in a position around xn, 
and for each x £ [0, 1], we calculate the limiting rescaled density function 
of &, < 6 < f, for the probability that the card in a position around 
xn will be a card with a number around bn. These density functions all 
have a discontinuity at x = b, and for each of them, the supremum of 
the density is obtained by approaching the discontinuity from one side, 
and, for certain values of the parameter, the infimum of the density is 
obtained by approaching the discontinuity from the other side. 



1. Introduction and Statement of Results 

Let Sn denote the symmetric group of permutations of [n] = {1, • • • ,n}. 
Our convention will be to view a permutation a G S'n as a deck of cards 
numbered from 1 to n and laid out in a row, where aj denotes the number 
of the card that is in the j-th position from the left. In a recent paper [4J, we 
analyzed the bias in the card cyclic to random insertion shuffle: remove and 

2000 Mathematics Subject Classification. 60C05, 05A05, 05A15. 

Key words and phrases, random shuffle, random permutation, total variation norm. 

1 



2 ROSS G. PINSKY 

then randomly reinsert each of the n cards exactly once, the removal and 
reinsertion being performed according to the original left to right order of the 
cards. The novelty in this nonstandard shuffle is that every card is removed 
and reinserted exactly once, unlike in any of the shuffles one encounters in 
the literature. The bias in this shuffle turned out to be surprisingly high, 
and possessed some interesting features. We describe one of these features 
now, the one that is the impetus for the present article. 

According to our convention, aj^ denotes the position occupied by card 
number j. Let Pn(id, •) denote the probability distribution for the above 
shufHe when it starts from the identity permutation. Let b G [0, 1] and let 
lim„_^oo bn = b, with 6„n being an integer in [n] . It was shown that the 
distribution function Fh{x) = lim„_>.ooPn(id, {cr^^^ < xn}) for the limiting 
rescaled position of a card with a number around bn possesses a density 
fb{x),0 < X < 1, which has a jump discontinuity; moreover, is constant 
and at its minimum to the left of the discontinuity, while as x approaches 
the point of discontinuity from the right, fbix) approaches its supremum. 
A similar phenomenon holds also for hx{b),0 < b < 1, the density for the 
limiting rescaled card number in a position around xn; more precisely, these 
densities also possess a discontinuity, the supremum is approached as b ap- 
proaches the discontinuity from the left, while for a certain range of the 
parameter value x, the infimum is approached as b approaches the point of 
discontinuity from the right. 

In the present paper we consider the above quantities for two cyclic ran- 
dom transposition shuffles. The first one works as follows: for j = 1, • • • ,n, 
on the j-th step transpose the card that was originally the j-th from the left 
with a random card (possibly itself). The method in this shuffle is simply 
the method of the above card cyclic to random insertion shuffle transferred 
from the context of insertions to the context of transpositions. Wc call this 
shuffle the card cyclic to random, transposition shuffle. The second shuffle 
works as follows: on the j-th step, transpose the card that is currently in 
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the j-th position from the left with a random card (possibly itself). We call 
this shuffle the position cyclic to random, transposition shuffle. 

These two shuffles are a lot more similar to one another than the above 
card cyclic to random insertion shuffle is to the corresponding position cyclic 
to random insertion shuffle, defined as follows: on step j, remove and ran- 
domly reinsert the card that is currently in the j-th position from the left. 
Indeed, it is easy to see that whereas by definition, every card gets removed 
and reinserted in the card cyclic to random insertion shuffle, in general many 
cards do not get removed and reinserted at all in the position cyclic to ran- 
dom insertion shuffle. On the other hand, it is easy to see that in both 
the card cyclic to random transposition shuffle and the position cyclic to 
random transposition shuffle, every card will get transposed. 

We will denote the probability measures on Sn corresponding respec- 
tively to the card cyclic to random transposition shuffle and the position 
cyclic to random transposition shuffle starting from a G by •) 
and p^°''((T, •)■ 

Theorem 1. i. Under p^'^ {id, ■), the random variable cr^^^, denoting the 
position of card number bnU, has the following behavior. A ssume that lim^_^Qo 
b G [0, 1] and that lim„^oo Xn = x & [0, 1]. Then 



fr\x) ^ hm npZ\^d, {a.\ = x„n}) = { 



e 



-1+6 



+ e-^ - e-^-^+^ X <b; 



ii. Under p'^^'^ {id,-), the random variable CTxnn, denoting the number of the 
card in position XnU, has the following behavior. Assume that lim^^oo = 
X G [0, 1] and that lim„_).oo bn = b & [0, 1] . Then 



= lim npX''\id,{a,„n = bnu}) 

n— >oo 



,-1+fe _^ ^-x _ Q-l-x+b ^ > 



X. 



Theorem 2. i. Under pn°^ {id, ■) , the random variable a f^^^, denoting the po- 
sition of card number bnU, has the following behavior. Assume that lim„-^oo &n = 
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b G [0, 1] and that lim. 



L„^.oo Xn = X G [0, 1] . Then 



fr{x)^ lim npr{id,{a, 



.-1 

bnU 



Xnn}) 



e 



+ e~^ X < 6; 



M. Under {id, ■) , the random variable (Jx^n, denoting the number of the 
card in position Xnn, has the following behavior. Assume that lim„_).oo Xn = 
X G [0, 1] and that lim„_>oo bn = b E [0, 1] . Then 



Remark 1. Of course, part (ii) of each to the theorems follows imediately 



from part (i), and gives h%^'^'^{b) = f^^'^'^{x) and hx°^{b) = /^°'^(x). Note also 
that it turns out that f^°\x) = h%^"^{b) and h^^b) = f§''"^{x). 



Remark 2. All four of the above densities have a discontinuity when the 
variable {x or b, depending on the density) is equal to the parameter (6 
or x). For the density f^^'^{x), the supremum is approached when the 
variable x approaches the value of the parameter b from the right. For the 
parameter 6 in the range [1— log 2, 1] pa [.307, 1], the infimum of the density is 
approached when the variable approaches the value of the parameter from 
the left. For the parameter not in this range, the density approaches a 
value less than 1 when the variable approaches the value of the parameter 
from the left, however the minimum of the density is attained at x = 1. 
For the density f^^^{x), the supremum is approached when the variable x 
approaches the value of the parameter b from the left. For the parameter b 
in the range [0, log2] [0, .693], the infimum of the density is approached 
when the variable approaches the value of the parameter from the right. For 
the parameter not in this range, the density approaches a value less than 
1 when the variable approaches the value of the parameter from the right, 
however the minimum of the density is attained at x = 0. Thus, under 
the card (position) cyclic to random transposition shuffle, the most likely 
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Figure 1. Card cyclic to random transposition shuffle: den- 
sity for limiting rescaled position of a card with a number 
around bn. 
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Figure 2. Position cyclic to random transposition shuffle: 
density for limiting rescaled position of a card with a number 
around bn. 

position for a card to end up in is a little to the right (left) of where it 
started, while the least likely position for it to end up in is either a little to 
the left (right) of where it started or at the very end (beginning) of the deck, 
depending on where it started. We conclude that the general phenomenon 
noted at the beginning of this paper with regard to the card cyclic to random 
insertion shuffle persists with the card cyclic and position cyclic to random 
transposition shuffles. See figures [l] and [2] respectively for graphs of f^^^'^ 
and f^"^- The graphs of h'^!^^'^ and hx"^ respectively are obtained from those 
of ff^°'' and f^'^'^'^ by switching the labels of the x and b axes. 

Let s and t denote generically b or x, and let fg{t) denote generically any 
of the above four density functions. The supremum of /* is e^^^'^ + e~^, 
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and the infimum is e 



-l+s 



+ e 



—s 



e . We have 



sup sup /* (t) = 1 + e 

0<s<l 0<t<l 



-1 



1.368, 



with the supremum approached as s and t both approach or 1, with one 
of the two variables strictly larger than the other one, the order depending 
on which of the four functions is considered. We have 



two variables strictly larger than the other one, the order depending on 
which of the four functions is considered. 

For the card cyclic to random transposition shuffle, let 



be respectively the expected limiting rescaled position for a card with a num- 
ber around bn and the expected limiting rescaled card number to be found 
in a position around xn, For the position cyclic to random transposition 
shuffle, let 



be respectively the expected limiting rescaled position for a card with a 
number around bn and the expected limiting rescaled card number to be 
found in a position around xn. 

Direct calculations give the following corollary. 

Corollary 1. i. 



The maximum of this function occurs at s = log 2 .693 with the value 
about .519, and the minimum of this function occurs at s = with value 
about .44^. 



inf inf fUt) = 2e-^ - e'^ f« .845, 



with the infimum approached as s and t both approach ^, with one of the 





e;Z\s) = EZais) = 1 - \e-'+^ - (1 + s)e-\ 
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II. 



1 



— s 



-1 




+ se 



The maximum of this function occurs at s = 1 with the value about .552, 
and the minimum of this function occurs at s = 1 — log 2 ~ .307 with value 
about .48I. 

In the case of the card cychc to random insertion shuffle, the density fb{x) 
degenerated when 6 — t- to include a 5-mass at x = of weight e~^. Using 
this as a starting point, we were able to show that the total variation distance 
between the distribution Pn(id, •) and the uniform distribution converges to 
1 as n — )• 00. Recall that the total variation norm between two probability 
measures and u on Sn is defined by 



In the cases at hand, we don't know whether the total variation distance 
between p^'^'^{id,-) and the uniform distribution or between pn°^{id,-) and 
the uniform distribution goes to 1. Let C/„ denote the uniform distribution 
on Sn- The above results just allow us to obtain the following extremely 
weak inequality. Let p*(id, •) be generic notation for either pn°^(id, •) or 



As was the case with the card cyclic to random insertion shuffle, the card 
cyclic to random transposition shuffle does not seem to have been studied 
before. On the other hand, the position cyclic to random transposition 
shuffle has been studied. From the point of view of mixing times, it was 
studied in P] and [3], where it was shown that as n — )• 00, the number of 
such shuffles needed to approach equilibrium is on the order logn. More in 
the spirit of this paper, the papers [S], [B] and [T] all studied various aspects of 
the distribution pn"^ (id, •), such as the limiting probability of a derangement 
occurring, or the limiting expected number of fixed points. In [5j it was 



H — vWtv = sup {fJ-{A) — iy{A)) 

AcS„ 




ires', 



'n 



^^^■''^(id,-). Then 
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shown that pn°^{id,a) > for all a G Sn, and that the inequality is an 

equality when a is the permutation which cycles every card to the right. 
In [T] it was also shown that for n > 18, but not for 3 < n < 17, the 
identity permutation has the highest probability; furthermore, as n — t- oo, 

p^°''(id, id) ~ — —fi ([5, p. 276], [Ij). By comparison, we note that 

in [3j we showed that for the card cyclic to random insertion shuffle, one has 



the sharp inequalities < p„,(id, ■) < where C„ = ~ 



is the n-th Catalan number. The results in [5], [6], and [1] do not at all allow 
one to determine whether or not lim„_5.oo ||Pn°'^(id, •) ~ ^«IItv = 1- It seems 
that the only result in those papers that could be used to give an estimate 
on the total variation norm is the one that says that under pn°^ (id, •), the 
probability of a derangement converges as n — t- oo to about .436 [6] , whereas 
under the uniform measure C/„ it converges to w .367, which gives an 
upper bound on the total variation that is even weaker than the one we 
obtained above. 

Question: Is it true that lim„_j.oo | |p$^^'^'^(id, ■)—Un\\TV = 1 and lim„_i.oo | |Pn°^(id, 

f/n||TV = l? 

Theorem [T] is proved in section 2 and Theorem [2] is proved in section 3. 

2. Proof of Theorem [T] 

Since p^^'^{id, {aj^ = a}) = pjf '^(id, {aa = j}), part (h) of the theorem 
follows immediately from part (i). We now prove part (i). 

For a,j G [n] with a ^ j, we consider p'^^'^{id, = a}), the probability 
that card number j ends up in position a. The shuffle has n steps, each of 
which is constituted by a transposition. One way for {cj^ = a} to occur 
is for card number j to move to position a on the j-th. step of the shuffle, 
and then for it never to move again. The probability of this occurring is 
-(1 — -)"^'' . The reader should convince himself that if on the other hand, 
card number j moves to position a on the j-th step of the shuffle, but is 
involved later on in another transposition, then it cannot end up in position 
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a. Another way for {crj^ = a,} to occur is if j < a, card number a is not 
moved before step a, on step a card number a is transposed with card number 
j, and then card number j is never transposed again. The probabihty of this 
is ^(1 — if i < ^- If 0^*3 the other hand, card number j is transposed 

again after step a, then it can not end up in position a. 

Assuming now that card number j is not moved to position a on the 
j-step or on the a-th step, we consider how else one can end up with the 
event {cr~^ = a}. The reader should convince himself of the following facts. 
(Remember that all the following statements are being made under the as- 
sumption that card number j does not move to position a on step j or on 
step a.) If card number a gets transposed before the a-th step, then card 
number j cannot end up in position a. (For example, say card number a gets 
transposed for the first time on the i-th step, with i < a. If card number i 
was not transposed before the i-th step, then on the i-th step, card number o 
was transposed with card number i. So now after the i-th step, card number 
i is in position a. On every later step k, card number k will be transposed 
with a random card. Since we are assuming that card number j was not 
moved to position a on the j-th step, there is no way for card number j to 
end up in position a. Similarly, if card number i was transposed before step 
i, say at step I, and card number I was not transposed before step I, then on 
the i-th step, card number a was transposed with card number I. And the 
same logic as above shows that card number j cannot end up in position a.) 

If card number a gets transposed for the first time at step a, and gets 
transposed with a card whose number is less than or equal to a, but not 
equal to j , then card number j cannot end up in position a. (The reasoning 
is similar to the above reasoning.) 

If card number a gets transposed for the first time at step a, and gets 
transposed with a card whose number I is greater than a, but then card 
number / gets transposed between step a and step I, then for reasons similar 
to the above, card number j cannot end up in position a. 
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However, if number card a gets transposed for the first time at step a, 
and gets transposed with a card whose number / is greater than a, and then 
card number / is not transposed between step a and step /, then there is 
stiU a chance for card j to end up in position a. One way would be if / > j 
and if on step /, card number / is transposed with card number j, and then 
after step I, card number j is never transposed again (this last requirement 
is possible because / > j)- Another way would be for card number / to get 
transposed on step I with card number m with m > I and m > j, then 
for card number m not to get transposed between step / and step m, then 
for card number m to be transposed with card number j on step m, and 
then for card number j not to be transposed again after step m (this last 
requirement is possible because m > j). Continuing to argue in this vein, 
we arrive at the following formula. 

Let Fk denote the first step on which card number k is transposed. Of 
course, from the definition of the shuffle one has <k. Let denote the 
last step strictly before the i-th step on which card number k was transposed. 
If there is no such step then define = oo. We will use the generic P when 
considering probabilities related to the random variables Fk and T^. Then 
we have 
(2.1) 

pr^(id,{a-i = a}) = -(1 - -r^^ + -(1 - -r-H,<,+ 



E E P{Fa = a)\J{P(T-=nAp(T^ 

m=l a=io<ii<---<im<n;im>j \l=l / 



■n+1 



We have P(F, = a) = (1 - P(7^; = ii.i) = ^(1 - ^Y'-'^-^-' and 

p(rpn+i ^ ^^-j ^ i^-^ _ l^n~im^ ^[^i^ last equality holding because 

im > j- Substituting this in (2.1) gives 

pT%d,{<Tj' = a}) = -(1 - -r-^ + -(1 - -)"-4,<a+ 
■' n n n n 

(2.2) n-a ^ ^ 

E\ ^ /j^ _\n—m—l/_\m+l 

^ n n 

m=l a=io<ii<---<im<n;im>j 
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If m > j—a, then the restriction im > j on the inner sum above is superfluous 
and we have 

E 

a=io<ii<--<im.<n;im>j 

li m < j — a, then we have 



n — a 
m 



E 1= E 

a=Jo<n<---<«m<n;«m>i «m=j + l 



ira-l - a 

m — 1 



Using this with (2.2) gives 



pr'(id,{cT-^ = a}) = ^ (1 - ^r-^ + +^ (1 - ^r-'ij<a+ 

■' n n n n 



(2.3) 



n—a 1 / 

EQ _ £^n-m-l/£^m+l (n — a 



m={j— a) ++l 



m 



+ 



Ef-y jj^yn— m— 1 /• jj^ ym+1 | \ ^ 
n n \ ^ 



m=l 



r=j+l 



r — 1 — a 
m — 1 



where any summation whose lower limit is greater than its upper limit is 
understood to vanish. 

Let ^Bin(7V,g) denote a binomial random variable with parameters N and 
q. We have 
(2.4) 



n—a ^ ^ 

\ ^ /j^ _\n—m—l/_\m+l 

^ n 



n 



n — a 
m 



m={j~a) + + 1 

Similarly, for r G {j + 1, • • • , n}, we have 



(1_ lp(X3 > (j-a)++l). 



(2.5) 



^ n n 

m=l 

S(l - r)"-^'^""'^(^Bin(.-l-a 1) < U - - !)• 



r — 1 — a 
m — 1 



n 



Now let a = Xnn and let j = bnn, with lim„_j>oo = x, lim„_j>oo = b, 



and b ^ X. We conclude from (2.4) that 
(2.6) 



lim n 

n— >oo 



EQ_£\n-m-l/£\m+l 1"^- XrJl 
n n \ m 

m={bnn—Xnn) + +1 



e-^(l-e-i+^), if X > 6; 
0, if X < b. 
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For all r G {hnn + !,••• , n}, we have 

(2-7) ^(^Bin(r-l-x„n,l) < (^n^ " ^n™)^ - 1) = 1, if X < 6. 



It then follows from (2.5) and (2.7) that 
(2.8) 

limn y (1- 



n 



, r=6„n+l 



r — 1 — a;„n^ 
m — 1 



" 1 1 

lim y -(l--)"-'^+^"'^-i = e'i-^' / e^dx = e"^ - e-^-^+^ if X < 6. 
n->oo ^-^ n n Jh 



If X > 6, then the left hand side of (2.8) is identically for large Ti- 



lt now follows from (2.3), (2.6) and (2.8) that 



(2.9) lim nprd(id, {a^^ = x„n}) = < 



e-i+^ + e"^, x>b; 



+ e 



□ 



3. Proof of Theorem [2] 



Since p^^{\d,{a~^ = a}) = pn°'^(id, {cJa = j}), part (ii) of the theorem 
follows immediately from part (i). We now prove part (i). 

The analysis here is similar to that in the proof of Theorem [T] so we will 
be less thorough here with the explanations. For a,j G [n] with a ^ j, 
we consider (id, {aj^^ = a}), the probability that card number j ends 
up in position a. The shuffle has n steps, each of which is constituted by 
a transposition. One way for = a} to occur is for card number j to 

be moved to position a on the a-th step and then for card number j never 
to move again. The probability of this is ^(1 — ^)"~"- Another way for 
{(jj^ = a} to occur is if a < j, card number j does not move before the j-th. 
step of the shuffle, on the j-th step it moves to position a, and then it never 
moves again after the j-th step. The probability of this is ^(1 — \)^~^ , if 
a < j. On the other hand, if card number j moves to position a on the 
a-th step and then moves again later on, it cannot end up in position a. 
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Similarly, if o < j, and card number j moves to position a on the j-th step 
of the shuffle, and then moves again later on, it cannot end up in position 
a. 

Assuming now that card number j is not moved to position a on the 
j-step or on the a-th step, we consider how else one can end up with the 
event = a}. By reasoning similar to that in the proof of Theorem jl| 

the only other way for this event to occur is if there exists an m > 1 and 
numbers satisfying j < ii < • • ■ < im 1^ n, with im > a, such that 

card number j does not move until step j, at step j it is moved to position 
ii, and then after that it moves again only at steps ii, I = 1, ■ ■ ■ ,m, with it 
moving at step ii, I = 1, ■■■ ,m — 1, to position and with it moving at 
step im to position o. The probability of the above occurring for a particular 
choice of m and {ii}]^^ is 

(1 — —y-'^ — h — iyi^J-i X • • • X —(1 — iv'n-^'^-i-iin _ }L\"--im _ 
n n n n n n n 

/-j^ jj^ Nn— m— 1 / jj^ -jm+l 

n n 

If m > a — j, then the restriction above that im > a is superfluous and we 
have 



n- J 
m 



E 1 

j<ii<---<im< >a 

If m < a — j, then we have 

n 

E 1= E 

j<ii<---<im<n;im>a im=a+l 

From the above analysis we conclude that 

pr(id, K' = «}) = -(!- -r-^ + -(1 - -r-'ia<j+ 

' n n n n 



im 1 j 

m — 1 



J 

n-j 



(3.1) 



^—^ n n \ m J 



m=(a— j)++l 

Ez-j^ _\n—m—lf_\m+l 
n n 

m=l 
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where any summation whose lower hmit is greater than its upper hmit is 



understood to vanish. Noting that the right hand side of (3.1) is the right 



hand side of (2.3) with the roles of a and j switched, part (i) of the theorem 
follows from part (i) of Theorem [T] □. 
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